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Abstract 

In the framework of the deformed quantum mechanics with minimal length, we consider the 
motion of a non-relativistic particle in a homogeneous external field. We find the integral rep- 
resentation for the physically acceptable wave function in the position representation. Using the 
method of steepest descent, we obtain the asymptotic expansions of the wave function at large 
positive and negative arguments. We then employ the leading asymptotic expressions to derive the 
WKB connection formula, which proceeds from classically forbidden region to classically allowed 
one through a turning point. By the WKB connection formula, we prove the Bohr-Sommerfeld 
quantization rule up to O We also show that, if the slope of the potential at a turning point 
is too steep, the WKB connection formula fall apart around the turning point. 
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I. INTRODUCTION 

One of the predictions shared by various quantum theories of gravity is the existence of 
an effective cutoff in the ultraviolet, i.e., a minimal observable lengthjl-g]. Some realizations 
of the minimal length from various scenarios are proposed. One of the most popular mod- 
els is the generalized uncertainty principle (GUP), derived from the modified fundamental 
commutation relation jl-l^] 

[X,P]=^(1 + /3P 2 ), (1) 

where f3 = (5 P p /h 2 = (3 /c 2 M 2 with the Planck mass M p = ^JhcjG and the Planck length 
£ p = a/ Gh/c 3 . /3q is a dimensionless parameter. With this generalization, one can easily 
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derive the generalized uncertainty principle (GUP) 



AXAP > + /3(AP) 2 }. 



(2) 



This in turn gives the minimal measurable length 



AX > A r 



%l p . 



(3) 



Eqn. ([T]) is the simplest model where only the minimal uncertainty in position is taken into 
account while the momentum can be infinite. 

In this paper we consider one dimensional quantum mechanics with the deformed com- 



mutation relation ([I]). To implement the deformed commutators ([T]), one defines 12|, Il5 | 



X — Xn, P — Pn 



V 

_ J n 



(4) 



where [Xq, Pq] = ih, the usual canonical operators. One can easily show that to the first 
order of /3, eqn. ([1]) is guaranteed. Henceforth, terms of O ((3 2 ) and higher are neglected in 
the remainder of the paper. For a quantum system described by 



H 



p2 

2^™ 



(5) 



the Hamiltonians can be written as 



H = H + H l + O (f3 2 ) 



(6) 



where H = £ + V (X ) and H l = ^-P 2 . Furthermore, one can adopt the momentum 
representation 

P, (7) 



Xq = ih—, P 

op 



or the position representation 



Xn — X, Pn 



h d 

i dx 



(8) 



The momentum representation is very handy in the discussions of certain problems, such 
as the harmonic oscillator 16|], the Coulomb potential 17|, |l8| and the gravitational well 19, 



20|. Recently, a wide class of problems, like scattering from a barrier or a particle in 
a square well 2lM23| are discussed in position representation. Moreover, in the position 
representation, it is much easier to derive and discuss WKB approximation in the deformed 
quantum mechanics analogously to in the ordinary quantum mechanics 2J] . Thus, we adopt 



the position representation in this paper. In the position representation, the deformed 
stationary Schrodinger equation is 

«*V (g) p dS>{x) 2m{E-V{x)) 

where we define Pa = ^h 2 (3 2 for later convenience. 

Although, the homogeneous field potential V (X) = FX is not studied so intensively as 
the quantum well, it has an important application in theoretical physics. In the ordinary 
quantum mechanics, the solutions to the Schrodinger equation with the linear potential 
are Airy functions, which are essential to derive the WKB connection formulas through a 
turning point. This motivates us to study the linear potential in the deformed quantum 
mechanics. 

In the deformed quantum mechanics with minimal length, the WKB approximation for- 
mulas are obtained in (2^]. In addition, the deformed Bohr-Sommerfeld quantization is used 
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261 ] . Therefore, it is 



to acquire energy spectra of bound states in various potentials 18, 
interesting to derive the WKB connection formulas through a turning point and rigorously 
verify the Bohr-Sommerfeld quantization rule claimed before, which are presented in our 
paper. Besides, we find that, if the slope of the potential is too steep at a turning point, the 
WKB connection algorithm fails around the turning point. This is not unexpected because, 
if one makes linear approximation to the potential around such a turning point for asymp- 
totic matching, the corrections to the wave functions due to the Hamiltonian Hi become 
dominant before one reaches the WKB valid region. 

This paper is organized as follows: In section [Til we give the integral representation of 
the physically acceptable wave function of the homogeneous field and its leading asymptotic 
behavior at large positive value of p. In section IIIIl we obtain the asymptotic expansions of 
the physically acceptable wave function at both large positive and large negative values of p. 
Section [TV] is devoted to deriving the WKB connection formula and the related discussions. 
In section [V] we offer a summary and conclusion. 

II. DEFORMED SCHRODINGER EQUATION 

Let us consider one- dimensional motion of a particle in a homogenous field, specifically 
in a field with the potential V (X) = FX. Here we take the direction of the force along the 
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axis of —x and let F be the force exerting on the particle in the field. As discussed in the 
introduction, the deformed Schrodinger equation for this scenario is 

In order to solve eqn. f flUj) . a new dimensionless variable p is introduced as 

p = (x - (2mF/h 2 ) * . (11) 

Eqn. (flOl) then becomes 

- aV 4) + ip" -pip = 0. (12) 

2 

where we define another dimensionless variable a 2 = £ 2 ^ (2mF/h 2 ) 3 and the derivatives are 
in terms of the new variable p. The linear differential equation (fT2|) is quartic and then 
there are four linearly independent solutions. We will shortly show that only one of them is 
physically acceptable. 



A. Physically Acceptable Solution 

The condition (3P 2 <C 1 validating our effective GUP model implies 

P (x\ P 2 |^> « (x \if,) =► a 2 W (p)\ « \if, (p)| • (13) 

This condition is also expected in the momentum space. Since the GUP model is only valid 
below the energy scale the momentum spectrum of the state \ip) should be greatly 

suppressed around the scale f3~^. It also leads to the condition ( 1T3"|) . Moreover, the condition 
(|T3l) and eqn. (fl2l) give 

a 2 |p|<l. (14) 

In other words, our GUP model, which is an effective model, is valid only when the condition 
( I14p holds. Considering that the Compton wavelength of a particle should be much larger 
than h\f/3 or lp in the GUP model, one can also obtain the condition (fl4|) in the classical 
allowed region where p < 0. In a field with the potential V(x), the kinematics energy of 
a non-relativistic particle is E — V (x) and its momentum is ^/2m (E — V (x)). Therefore, 
the fact that the Compton wavelength of the particle A c = -, h = is much larger than 

y/2~m(E-V(x)) 

£/3 yields a 2 \p\ 1. In the remainder of our paper except subsubsection HV B 1\ we assume 
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a < 1 which is useful to derive WKB connection formula around a smooth tuning point. 
One needs to consider a ~ 1 scenario only when it comes to the WKB connection around a 
sharp turning point. 

We notice that E < V for p > 0. The wave function ip is then exponentially damped 
for large positive value of p. Thus, one needs to evaluate asymptotic values of ip (p) at large 
positive value of p to find physically acceptable solution to eqn. ( Fl2|) . Note that, only when 
a <C 1, one can analyze asymptotic behavior of ip (p) at large positive value of p in the 
physically acceptable region where a 2 \p\ <C 1. To determine the leading behavior of ip (p) 
at large positive value of p, we make the exponential substitution ip (p) = e s ^ and then 
obtain for eqn. (ITS]) 

s" + s' 2 -p-a 2 [s^ + 6s' 2 s" + 3s" 2 + 4sV 3 > + s' 4 ] = 0. (15) 



Eqn. (fl5l) is as difficult to solve as eqn. (112]) . Here our strategy to find the asymptotic 
behavior of %fj (p) from eqn. ( Tl5l) is as follows [2 71]: 



(a) We neglect all terms appearing small and approximate the exact differential equation 

with the asymptotic one. 

(b) We solve the resulting equation and check that the solution is consistent with approxi- 

mations made in step (a). 

It is usually true that higher derivative terms than s' are discarded in step (a). Therefore, 
we reduce eqn. (fl5|) to the asymptotic differential equation 

s '2_ a 2</4^ (16) 

Solving eqn. (fl6|) gives four solutions for s', two of which are discarded considering (3P 2 1. 
Taking asymptotic relation ffT3l) into account, one can further reduce the quartic equation 
f|T6l) to a quadratic equation 

s' 2 ~ p. (17) 
which has only two solutions for s'. The two solutions are s' ~ i^/p, and, therefore, 

ip (p) ~ exp ^±-p2^ , at large positive value of p, (18) 
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where — is for the physically acceptable solution. It is easy to check that the solution 
s' ~ ±-y/p satisfy the assumptions 

s",s' 2 s", s^,s" 2 , s's™ and s< 4 > « p, 

as long as p 1. 

It is interesting to note that the two discarded solutions of eqn. ( !T6|) are 



s> ~ ±J£ ^= , (19) 

\/2a 



which become s' ~ ^775 when a 2 p <C 1. The resulting wave functions are ijj (p) ~ 
exp N=^=Y They are not physical states since they fail to satisfy the condition (Tl3]) . One 
can also see that these two solutions are discarded according to the low-momentum consis- 
tency condition in 28]. In summary, assuming a <C 1, we find that the leading asymptotic 
behavior of the physically acceptable solution is exp ^— |p^j for p ^> 1. In addition, we 
only analyze the solution in the region a 2 \p\ <C 1 where the GUP model is valid. 

B. Integral Representation 

The differential equation eqn . ffl2l can be solved by Laplace's method. Please refer to 



mathematical appendices of [29j for more details. Define the polynomials 

P(t) = -aH A + t 2 , Q(t) = -1, (20) 

and the function 

Integral representations of the solutions to eqn. ffT2|) are then given by 

$ (p) = - / exp (pt) Z (t) dt 
Jc 

a 2 t 5 t 3 



exp (pt + — - - dt, (22) 

where the contour C is chosen so that the integral is finite and non-zero and the function 

( a 2 t 5 t 3 \ , . 

V(t)=exp[xt+ — --), (23) 



vanishes at endpoints of C since the integrand of eqn. (1221) is entire on the complex plane 
of t. Now that exp (xt + — ~ exp ^^-^-jfor large t, we need to begin and end the 
contour C in sectors for which cos 56 < (setting t = \t\e ie ). There are five such sectors, 
specifically 

9eB n = g— ^, g— ^- , n= 0,1,2,3,4,5. (24) 

Therefore, any contour which originates at one of them and terminates at another yields a 

solution to eqn. (fl2"|) . One could then find four linearly independent functions of the form 

r ( a 2 t 5 t 3 \ 
h{p)= / exp [pt + — - - dt. (25) 



The asymptotic expression for Jj (p) for large values of p is obtained by evaluating the 
integral eqn. ( |25i) by the method of steepest descents. 

III. ASYMPTOTIC EXPANSION 

First we briefly review the method of steepest descent to introduce some useful formulas. 
This technique is very powerful to calculate integrals of the form 

I(p)= [ g(z)ePMdz, (26) 



Jc 

where C is a contour in the complex plane and g (z) and / (z) are analytic functions. The 
parameter p is real and we are usually interested in the behaviors of / (p) as p — > ±oo. 
The key step of the method of steepest descent is applying Cauchy's theorem to deform 
the contours C to the contours consisting of steepest descent paths and other paths joining 
endpoints of two different steepest descent paths if necessary. Usually, the joining paths are 
chosen to make negligible contributions to / (p). It is easy to show that Im/ (z) is constant 
along steepest descent paths. When a steepest descent contour passes through a saddle point 
zq where /' (20) — 0, / (z) and g (z) are expanded around zq and Watson's lemma is used to 
determine asymptotic behaviors of / (p). Specifically, consider a contour C through a saddle 
point Zq. A new variable r is introduced as r = f (z) — f (zo) to calculate I (p). The saddle 
point Zq divides the contour C into two contours C\ and C2. Generally, r monotonically 
increases from —00 to zero along one contour, say G\ and monotonically decreases from zero 
to —00 along C 2 . Thus, the integral becomes 



I (p) = exp [pf (z )} 



f° dz f^°° dz 

J g(r)exp[pT]—\ Cl dT + J g (t) exp [pr] —\c 2 dr 



(27) 



The physically acceptable solution can be represented by an integral 

a 2 t 5 t 3 ~ 



Hp) 



J exp ^pt + 



dt, 



(28) 
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where C is any contour which ranges from t = exp ( — ^p) oo to t = exp (^p) oo. In fact, 

as we show later in the section for positive p, there exists a steepest descent contour from 

t = exp (— 2p) oo to £ = exp ( 3 J 2 ) oo, which C can be deformed to. Moreover, the integral 

on such a steepest descent contour yields the required asymptotic behavior of I (p) at large 

positive value of p. Here the exponent in the integrand has movable saddle points. Making 

i 

the change of variables s = \p\ 2 1, one gets 



I(P) = \P\ 



exp m oo 



cxp(-22l loo 



exp 



as 5 s 3 

±s H 

5 3 



ds 



\P 



2 / exp 

Jcxn(-: 



\p\ 2 f±(s) ds, 



'exp(-2ji)c 

where + for p > and — for p < and a = a 2 \p\ <C 1 in the physical region. 



(29) 



A. Large Positive p 

For p > 0, we have 



as 5 s 3 



/+(*) = * + — -y (30) 



There are four saddle points given by f' + (s) = at 

Vl-4a + Vl-4a 

s = ±A + = ± = and s = ±rj + = ± = . (31) 

V2a V2a 



Our goal now is to find a steepest descent contour emerging from s = exp (— ^p) oo to 
s = exp (^p) oo. We will show that such a contour passes through s = — A+. To find the 
contour we substitute s = u + iv and identify the real and imaginary parts of /+ (s) 

1- — + — + v 2 - 2cm V + av 
3 5 

+iv ^1 - m 2 + aw 4 + y - 2auV + . (32) 
Since Im/ + (— A + ) = 0, the constant-pahse contours passing through s = — A + must satisfy 

v 1 - u 2 + au 4 + — - 2au 2 v 2 + — = 0. (33) 



Therefore, one of the constant-phase contours passing through s = — A + is 

C : != \ 1 + lav 2 — \ 1 — 4a + -av 2 H — -a 2 v 4 + iv , for — oo < v < oo, 

2a V V 3 5 



which is a steepest descent contour. In fact, around the saddle point s = — A + , one finds on 
the contour C 

s ~ -A + + bv 2 + iv, (34) 

and hence, 

/+ (s) = f + (-A+) - V ^fl (-A+) + (t; 3 ) , (35) 

where 6 is a positive real number. Since /" (— A+) is real and positive, the contour C is 
indeed a steepest descent contour. Note that C goes to s = exp (— 2p) oo as f — >■ — oo and 
s = exp (^p) oo as f — 7- oo. In order to evaluate asymptotic expansion of I (p), we break up 
the contour C into C\ and C2, corresponding to above and below of s = — A+. Define 



r = /+ (a) - / + (-Ah 



(36) 



where r monotonically decreases from zero to —00 as one moves away from s = — A + along 
Ci to s — exp (^p) 00 and along C 2 to s = exp (— ^p) 00, respectively. Since f + (— A + ) = 0, 
the expression for s in terms of r can be expressed as a power series of y/—r. Then, noting 
that — r = (±V — r ) 2 ; one nas 



-A + + ^a,(±v^) : 
i=i 



(37) 



where a^ can be obtained by substituting eqn. (J3"7]l into eqn. fl36|) and equating powers of 
V'— r on both sides of the equations. It is easy to find 



d\ = Ia 



(38) 



where one finds Ima! > 0. The contour C\ is in the second quadrant and hence, + sign is 
chosen in eqn. fl37|) for C\. Therefore, 

ds 



I f 


5 t 




p 2 / exp 


P 2 f 


As) 


JCi 







exp (p 2r ) 3~^ r 



exp 



' 2 / + (-a + ) x; 



/ j 3j-2 



^ r(l). 



(39) 
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For the contour segment C2, the sign of yfr occurring in eqn. ( 137)) has to be reversed. 
Moreover, the limit of integration on C2 in the variable r ranges from —00 to 0. Thus, 



1 

P 2 



/ exp 


3 „ 
P 2 f 


4-00 


ds ~ — exp 


P 2 /- 


h(-A+)] 















t 2p^ V2 



r ■ (40) 



Combining eqn. (139)) and eqn. ( 140|) . we easily find 



/ (1 « p « O - ^ r l "" JJ E ( 3 + i) . (4i) 



P 3 i=o P 2 

B. Large Negative p 

As for p < 0, the exponent in the integrand of / (p) is 

f-(s) = s + ^-j. (42) 
Thus, one as well finds four saddle points given by f' + (s) = 

. = ±A_ S ±^ 1 -^ + 4a and . = ±„_ S ± Vl + ^ +4 °. (43) 



As before, our objective is to find steepest descent contours passing through the saddle 
point(s) in eqn. ()4"3)) that emerges from s = exp (— ^p) 00 to s = exp (^p) 00. Substituting 
s = u -(- if, we obtain the real and imaginary parts of /_ (s) 



(w 2 CLU^ 
-1 - — + — + w 2 - 2cmV + at)' 

+iv ^-1 - w 2 + cm 4 + y - 2au 2 v 2 + . (44) 

We have already shown that only one steepest descent contour passing through s = — A + is 
sufficient to evaluate asymptotic behavior of / (p) for large and positive p. However for large 
and negative p, things are a little bit more complicated. Instead of one steepest descent 
contour, it turns out that we need three steepest descent contours passing through ±A_ and 
r]_, respectively, to connect two endpoints at s — exp (±^p) 00. 

First consider the steepest descent contour through s = — A_. Since / + (— A_) is a pure 
imaginary number, the steepest descent contour must satisfy 

(V 2 (XV 4 \ 

-1 - u 2 + au 4 + — - 2au 2 v 2 + — ) = f- (-A-) . (45) 

11 



Solutions to the last equation give us a constant phase contour C-\_ passing through s = 
— A_, which emanates from s = exp (— ^p) oo and finally approaches s = exp (— y) oo. The 
contour C_a_ actually is composed of three segments as 

1 



A-,l 



A_,2 



1 + 2av 2 



F-X- (v) + iv, for — oo < v < — Im A_ 



x 2a 



1 + 2av 2 — \ F_\ (v) + iv, for — Im A_ < v < —v , 



where we define 



1 + 2av 2 + \ F-\ (v) + iv, for — t> > v > — oo, 



16 



F± x _ (v) = 1 + 4a + ^av a + — -aV + 



4a/_ (±A. 



3 5 iv 

and f is a solution to F + \_ (v) = that satisfies < % < 1. It is straightforward to 

verify that, along C_a_, Re /_ (s) monotonically increases from — oo to as one moves from 

s = exp (— ^p) oo to s = — A_ and then monotonically decreases from to — oo as one moves 

away from s = — A_ to s = exp (— y) oo. Hence, the contour C_a_ is indeed a the steepest 

descent contour passing through s = — A_. Now we calculate the contour integral on C-\_. 

Introduce 

r = /_(«)-/_ (-A_) , (46) 
which r is real on C_a_ an d varies from — oo to zero and then to — oo along C-\_. Then, 
one has 

oo 

S = -A_ + $>(±v^) J , (47) 
i=i 

where 6j can be obtained by substituting eqn. (j4"6"|) into eqn. (]4"T]). One easily gets 

'7T 



b\ = exp ( - 



(48) 

Since Re exp > 0, one has —\J—t for C_a_,i and a/— r for C_a_,2 + C_a_,3 in eqn. fj47|) . 
Therefore, 







p 2 / exp 


H 1 /-(*)] 


</c_ A _ 





2p2 exp 



exp |p| 2 ,/L(-A 



exp (|p| 2 rj J2 (2J + X ) ^ (v 7 ^) ^v 7 ^ 
i=o 



(49) 



(2j + l)6 



IP 



i=o IP 



12 



Analogously, one can readily write down a constant phase contour C + \_ passing through 
s = — A_, which starts from s = exp (y) oo and ends at s — exp (^p) oo. As before, C + \_ 
consists of three segments 



a 



+A-.1 



1 + 2at> 2 + -v / F + a_ (f ) + iv , for + oo > v > v , 



c. 



1 



A_,2 



a 



2a 
1 



1 + 2av 2 — \ F +x _ (v) + iv, for v > v > Im A_ 



+A_,3 



1 + lav 2 — \ F +x _ (v) + iv, for Im A_ < v < +oo. 



It is also straightforward to verify that C + \_ is a steepest descent contour as well. Setting 

r = /_ («) - /_ (A_) , (50) 

one finds r is real on C + \_ and varies from — oo to zero and then to — oo along C + \_. Note 
that /+ (s) is an odd function and X*_ = — A_. Taking complex conjugate of both sides of 
eqn. ( |50l) . one then has on C+\_ 



3=1 



(51) 



Since Re b\ > 0, one has \f—r for C + a_,i + C*+a_,2 and —\J—t for C + a_,3 in eqn. (I5T]) . 
Therefore, 







ds ~ 


p 2 / exp 


)p\h-(s) 


Jc +X _ 







exp 



|p| 5 /-(A- 



\p 



E a, r( J + -). (52) 



Since the values of Im/_ (s) are different on C±\_, it is obvious that we need a third 
contour which joins C±\_ up at s = exp (±y) oo, respectively. Here, we consider a constant 
phase contour C r; _ connecting s = exp (— oo to exp (y) oo that passes through Since 
Im / (s) = Im/_ (r^_) = on the contour C rj _, one finds 



1 



8 16 

1 + 2av 2 + \jl + Aa H — at> 2 H a 2 v A + iv, for — oo < v < +oo, 



x 2a V V 3 5 

is a curve of steepest descent. On C v _, define 



T = f- (S) ~ f- (V-) 



(53) 
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which is real on C v _ and varies from — oo to zero and then to — oo along C v _. Then, one 
finds 

oo 

s = \_ + J2c j {±V=f) j , (54) 
i=i 

and 

ci = ! vS5 (55) 

Similarly, we break up the contour C v _ into C V t i and C V t 2, corresponding to above and 
below of s = 77_ with a/— t for G n _^_ and —y/^r for C n _^ in eqn. ( |54l) . Thus, 









p 2 / exp 


)p\h-(s) 


ds ~ 


J Cr\_ 







exp 



\p\~ 2 f-(v-) 



IP 



E 



|p| 2 V " 



(56) 



Note that although paths C v _ and C±\_ never join up at s = exp (±y) oo, the integrand 
exp [/_ (s)] ~ exp ^^y- s 5 ^ tends to zero exponentially. Therefore, there is no contribution 
from a connecting path from C v _ and C±\_ at a distance R from the origin in the limit 
R — > oo. As a result, the integral / (p) equals to the sum of three contour integrals on the 
different steepest descent curves C v _ and C±\_. Combining eqn. (T49|) . eqn. f )52|) and eqn. 
(|56|) gives the full asymptotic expansion of / (p) for large and negative p 



I (-1 > p > 



-a 



00 Im ( exp 

! )~2<£- 

i=o 



-HV-(A-) 



^fii+i) r(j . 



+ 



\p\ 4 Ipl 2 

exp ~ (2j + l)c 2i+1 



i=o 



\P\ 2 



(57) 



IV. WKB APPROXIMATION 



The authors of 24J find the WKB approximation in deformed space with minimal length. 



In 



24J, they consider the deformed commutation relation 

[x,p]=ihf (P), 



(58) 



where / (P) is an arbitrary function of P. In our paper, we set / (P) = 1 + (5P 2 . Defining 
Pip) 



dP (p) 
dp 



f(P), 



(59) 
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and p (P) an inverse function of P (p) , they find the physical-optics approximation to the 
solution of the deformed Schrodinger equation 



nl 'h d\ . . . 2m(E-V(x)) . , s 



0. 



is 



ip (x) 



1 



Ci exp 



pdx 



+ C 2 exp 



pdx 



(60) 



(61) 



VpfiP) 

where P = ^2m (E -V(x)) in eqn. (16T1) . It is also shown there that, if eqn. (I6TT) is valid, 
the condition 



p 2 > h 



dx 



(62) 



has to be satisfied. However, the condition eqn. ( 162]) fails near a turning point where 
P = 0. Thus, if we want to determine bound state energies, we need to be able to match 
wave functions at the turning points. Here we considers a potential V (x) with its classical 
turning point located at x — 0. A linear approximation to the potential V (x) near the 
turning point x = is 

V (x) w V (0) + Fx, (63) 

where F — V (0). The linearized potential (1631) is discussed in the previous two sections. 
Our discussion shows that the parameter a = £p (2m \F\ /ft 2 ) 5 plays an important role in 
analyzing asymptotic behaviors of the solutions. When a< 1, the physically acceptable 
solution can exist at large argument p while the condition ( I14p still holds. Accordingly, a 
turning points is called a smooth one if a = £p (2m \F\ /fr 2 )* <C 1. Otherwise, it is called a 
sharp turning point. 



A. WKB Connection through a Smooth Turning Point 

Now we want to match WKB wave functions at a smooth turning point in the deformed 
space with / (P) = 1 + (3P 2 up to O ((3). Suppose x = is a smooth turning point, which 
means V (0) = E, and V > E for all x > 0. The region to the left of the turning point is 
classically forbidden where the wave function must be damped and become zero at infinity. 
Thus, far from x = 0, the wave function has the form 



ip (x) 



1 



VpTTp) 



C exp 



pdx 



for x > 0. 



(64) 
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To the right of the turning point, the wave function is given by 
1 



if) (x) 



C\ exp 



pdx 



+ C 2 exp 



h 



pdx 



for x < 0. 



(65) 



Around the turning point, x is small and P ~ \/2mF \/ —x. In this region, we may approx- 
imate eqn. ( 1641) and eqn. ( 1651) by 



2 x§ ' 1 + Yo a + O ( a2 ) 



3h 



for x > 0, 
(66) 



if) (x) 



C\ exp 



f C 2 exp 



| W t(l-Aa + (a*) 



for x < 0. 
(67) 



The criteria (162]) for validity of the WKB approximation is satisfied if 



(68) 



where we neglect f3P 2 in derivation. On the other hand, when the potential is linearized 
around the turning point x — 0, the Schrodinger equation becomes 

d 2 ip (x) 2 d 4 ip (x) 1m\xx 



dx 2 



p dx 4 



h 2 



-iP (x) « 0, 



(69) 



34 



where /3 = t>t§. To solve the approximate differential equation, we make the substitution 



p = x (2mF/h 2 )' 



(70) 



In terms of p, the solution to eqn. (159*1) which matches eqn. fl66|) and eqn. (1671) in two 
different limits is actually I (p) calculated in the section IIHI Specifically, the solution is 



if; (x) = DI (p) = DI (x (2mF/H 2 ) , 



(71) 



where D is a constant to be determined by asymptotic matching. It is easily shown from 
( |68|) that there exists overlap regions where both WKB approximation and eqn. (169]) hold. 
In the overlap regions, one finds |p| 3> 1 and \x\ <C 1. Therefore, we approximate I (p) by its 
leading asymptotic behaviors for large argument in the the overlap regions. The appropriate 
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formulas are 



I (l< p < a 2 ) 



/ (-1 > p > -a" 2 ) 



cxp 



2pt 



2i0F(l- f + 0(a 2 )) 



sin 



2 |p|t 



3a 
10 

3a 
10 



+ O (a 2 ) 



0(a 2 ) 



(72) 



(73) 



where a < 1 for a smooth turning point and a = t 2 g {2mF / h 2 ) z \p\ C 1 as required by the 
condition f[T5]) . Requiring that eqn. ( 172"]) and eqn. (175]) match eqn. (155]) and eqn. f ]57]) 
in the overlap region, respectively, gives Ci = — iCe m / A and C*2 = iCe l7r ^ 4 up to 0(/3). In 
summary, in the overlap region, we find WKB solutions and the asymptotic values of the 
solution to the Schrodinger equation with a linear approximation to the potential V(x). 
Then, by making eqn. (172]) and eqn. (175]) match eqn. ( l66j) and eqn. (|67j) respectively, the 
WKB connection formula with the deformed commutator [X, P] = %h (1 + /3-P 2 ) is obtained 
up to O ((3). The connection formula around a smooth turning point is put in a way that 



C 



cxp 



pdx 



2C 



sm 



7T 



prfx + - , up to O {13) , (74) 



which is directional, just as in ordinary quantum mechanics 27j] . The analysis always 
proceeds from classically forbidden region to classically allowed one. For bound states, 
the uniqueness of the wave function in the classically allowed region leads to the Bohr- 
Sommerfeld quantization condition 

nb 



pdx 



i) H — ] nh, up to O ((3) , 



(75) 



where a and b are two smooth turning points for the potential V (x). Notice that although 
eqn. ( 175]) is claimed in [24j, one still needs to obtain the connection formula to derive eqn. 
r5]) rigorously, which is not presented in 24]. 



B. Discussion 



1. Sharp Turning Point 



Near a sharp turning point x = 0, not only the WKB approximation falls apart but also 
matching the two WKB solutions across the turning point stops making sense. In fact, from 



17 



the previous subsection, one finds that the asymptotic matching is valid as long as there 
exists an overlap region where 1 <C \p\ <C a~ 2 . However, such region doesn't exist unless 
a < 1, which means that the asymptotic matching fails through a sharp turning point. 

It can be shown, through (|62p . WKB approximations are valid as long as \x\ ^> 
(2m |F| /h 2 )~* in the region where the potential is approximated by a linear one. Put 
another way, if there exists a region where both WKB and linear approximations are valid, 
one finds \x\ ^> (2m \F\ /h 2 ) 1 for such a region. When \x\ ^> (2m \F\ /h 2 ) 5 , we have 

\PP 2 \^ ^ t r»l (76) 

(2m\F\/h 2 )* (2m\F\/h 2 )* 

for a sharp turning point. However, \(3P 2 \ <C 1 is required by the GUP model. This means 
that, as moving away from the sharp turning point, one is far beyond the region where the 
linear approximation to the potential is good before reaching the WKB valid region. One 
might resort to a higher order approximation to the potential and asymptotic matching in 
the overlap region to find WKB connection formula through a sharp turning point. 



2. O ((3) vs. O (ft) 

When h can be regarded as a small quantity, the approximate solution to the deformed 
Schrodinger equation 



d 2 ip(x) 2h 2 (5d A i)(x) 2m (E — V (x)) , . 

H W {x) = 0, 



(77) 



dx 2 3 dx 4 h 2 

is easy to find using WKB analysis. To be specific, the approximate solution is expressed in 
an exponential power series of the form 



if) (x) ~ exp 



^ oo 

- K n Sn (X) 



n=0 



The authors of l24l finds 



S 1 = --\n\2Pf(P)\. 



Since here / (P) = 1 + (3P 2 , we have for Si 



Si w In 







\P\ 



p 2 + o (p 2 ) 



(78) 



(79) 



(80) 
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Moreover, the leading order (in terms of 0) of the S 2 is just the WKB O (h 2 ) correction 



calculated in the ordinary quantum mechanics. Therefore, we obtain 



29] 



S ^&* + J ^ dx+ (81) 
If one uses WKB approximations to evaluate quantum gravity induced corrections, say to 
energy levels or tunnelling rates, one may want to have 

PP 2 > hS 2 . (82) 

Otherwise, the quantum gravity correction (~ O (/3)) on the first order WKB approxima- 
tion (~ O (h )) could be overwhelmed by the second order WKB approximation (~ O (h)). 
Suppose a is the characteristic length of the potential V (x) , for example the width of a 
square-well potential. Then we can get a rough estimate on S 2 

hS 2 ~\~ -, (83) 
aP a 

where A is the de Broglie wavelength of a particle with momentum P. As a result, the 
condition ( I8"2"j) becomes 

It is interesting to note that the condition (151]) is a rough estimate and a more accurate 
estimate could be obtained once the form of the potential is given. 

Taking into account the constraints ( BID on the de Broglie wavelength A of a particle, 
one may conclude that the WKB approximation is not a powerful tool to calculate quantum 
gravity corrections unless the energy of the particle considered is high enough. However, 
there is an exception if the corresponding Schrodinger equation in the ordinary quantum 
mechanics can be solved exactly. In this case, O (/3) corrections calculated on the WKB first 
order approximation are just quantum gravity corrections to exact results up to O (/3) O (h°) 
even without having (184p required. For example, if we employ WKB analysis to calculate the 
energy spectrum of a bound state in the deformed space, the energy levels can be represented 
by a series in powers of H 

E n = J2VE nJ (P) } (85) 

j=0 

where E n j ((3) can be expanded in terms of j3 

00 

E n AP) = Y,P kE n,r ( 86 ) 



k=0 



19 



If on the first order WKB approximation, one calculates E n $ (j5) up to O (f3) 

E nfl 09) = E% + (3El fi + O (/3 2 ) , (87) 

the energy levels are 

E n = E% + (5E l nfi + O (/3 2 ) +0(h). (88) 

In order to have eqn. (188]) make sense, one requires PE^ > 0(h). On the other hand, if 
we know the exact result E n with = 0, namely -E 1 ^ 

oo 

i#»> = ( 89 ) 
i=o 

eqn. (1551) becomes 

£ n = + /3< + O (ft) (/3) + O (/3 2 ) . (90) 

Since O (ft) O (ft) is automatically smaller than /3Eq, eqn. (l9T)j) always makes sense as long 
as 0(h) < 1. 

To illustrate our points, we use the WKB approximation to derive the energy levels of a 
particle confined to the one-dimensional potential V (x) — F \x\ whose turning points are 

E , E 

a = --,b = -. (91) 
The energy quantization condition ( 17"5|) from first order WKB approximation then becomes 

f - % /_* v I - Mdx - S- t (f _ w ) ' dx = v + v * + ° m • (92) 

where lp — (h 2 /2mF) 1 is the characteristic length of the potential V (x) = F \x\. From the 
last equation, we obtain 



F 



(l + ^ + 0(/3 2 )+O(ft)), (93) 



z 

where £ n = £p [§ (n + |) 7r] 3 . What is (9 (ft)? The second order generalization of eqn. ( 1751) 



with /3 = is given in [27] 

3 /' 
2 / 



5 / 1 v ^ - NdI + S = (" + * + ° (R2) • (94) 



which gives 

Ek 0) . ( . a 



(i - 4 + ° (fi2) ) ■ <95) 
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We can then estimate O (h) through ([95 



0(h)~% (96) 

which can also be easily obtained by dimensional analysis. If one wants the first order 
approximation fl93l) to make sense, the second term in ([931) should be comparable to or 
larger than O (h) and then one gets 

(97) 




The de Broglie wavelength of a particle with energy E n ~ F£ n is 



h 



X n ; = jL. (98) 

^FT n 
Thus, the inequality ([57)1 reads 

K<l F {£)\ (99) 

which is much milder than (I84p . In a practical way, h and can be expressed in terms of 
£p, If and l n . In fact, it is easily shown that 



V. CONCLUSIONS 

In this paper, we considered a homogeneous field in the deformed quantum mechanics 
with minimal length. The physical motivation for this is to obtain the WKB connection 
formula and prove the Bohr-Sommerfeld quantization rule rigorously in the deformed quan- 
tum mechanics. By studying the leading asymptotic behavior of the physically acceptable 
wave function in the physical region, we found the contour for its integral representation. 
Through the integral representation, the asymptotic expansions of the physically acceptable 
wave function at both large positive and large negative values of p were given. 

Then, we used the obtained asymptotic expansions to get the WKB connection formula, 
which proceeds from classically forbidden region to classically allowed one through a smooth 
turning point, and had the Bohr-Sommerfeld quantization rule proved rigorously up to 
O (f3). A new interesting feature appearing in the presence of deformation is that our WKB 



21 



connection formula do not work for a sharp turning point. The connection through such a 
point might need a higher order approximation to the potential near it. 

Finally, we discussed the competition between the quantum gravity correction on the 
first order WKB approximation and the second order WKB approximation. If the former is 
not overwhelmed by the latter, the energy of the particle considered should be high enough 
according to (I84p . We also showed that, if the energy levels of a bound state are given 
in the ordinary quantum mechanics, the deformed energy levels are 

E n = + pE x nfi + 0(h)0 + O (/3 2 ) , (101) 

where P>E\ is the O ((3) quantum gravity correction on the first order WKB approximation. 
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